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Thermodynamic Properties of Quantum 
Mechanical Gravitational plus Harmonic Oscillator 
Potential using the Proper Quantization Rule 


Benedict | Ita, Nelson A Nzeata, Hitler Louis 


We obtained by using the proper quantization rule, the analytical solutions of the Schrédinger equation for the Deformed Harmonic 
Oscillator potential. The energy levels of all bound states as well as the normalized wave-functions are obtained. The vibrational 
entropy, heat capacity, partition function, mean free energy and, chemical potential were analyzed to further investigate the behavior 
of the mixed potential under investigation. 


INTRODUCTION 

The exact solution of both relativistic and non-relativistic equation for most of the physical potential possesses difficulties in area where the 
application of the approximation becomes necessary. i.e, in the case of the Schrédinger equation when the arbitrary angular momentum quantum 
number I is present, one can only solve the non-relativistic equation approximately using a suitable approximation scheme [1]. Some of such 
approximations which has successfully produced results include conventional approximation scheme proposed by Greene and Aldrich [2], improved 
approximation scheme by Jia et al. [3], elegant approximation scheme [4], the pekeris approximation which is based on a Taylor expansion series of 
the centrifugal potential [5] and recently a good approximation by Yazarloo et al. [6]. Some of these potentials include the Manning-Rosen potential 
[7], Rosen-Morse potential [8], Mobius square potential [9], Shifted Deng-Fan potential [10], Shioberg potential [11], Morse and generalised Morse 
potential [12].The Manning—Rosen potential have been widely investigated both in its pure state and when mixed for relativistic and non-relativistic 
equations using various methods where different approximation scheme were employed in the literature [13-23]. 

In this present work, a mixed potential known as the Quantum Mechanical Gravitational plus Harmonic Oscillator Potential, is investigated 
through the proper quantization rule (PQR). Recently, the thermodynamic properties of different potentials have been studied. Oyewumi et al. [24] 
studied the thermodynamic properties of the Schrodinger equation with a shifted Deng-Fanpotential via AIM method. Larkin and Filinov [25] 
investigated the thermodynamics of relativistic Newton—Wigner particle in external potential field. Suparmi et al. [26] studied the thermodynamics 
properties of diatomic molecules with q-deformed modified Poschl-Teller plus Manning Rosen non-central potential in D dimensions using 
SUSYQM approach. Sincethen, different authors have investigated the thermodynamic properties for some physical systems [27, 28] 

We organize our work as follows; firstly, we give a brief introduction in section one. Secondly, in section two we discussed briefly the proper 
quantization rule. In section three, we obtained the approximate solution of the Schrodinger equation for any |-state. Section four will be dedicated to 
the thermodynamic properties of the MRHP model. Section four will present results and discuss the work. Finally the conclusion will be given in the 
last section. 


PROPER QUANTIZATION RULE 
The proper quantization rule has been used to investigate different quantum potential to obtain both exact and approximate solutions of non- 
relativistic particle equation by Qiang—Dong [29]. For the non-relativistic equation, the one dimensional Schrodinger equation expressed as 


0) = SRE - UO) () 
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Can be written as 
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where @(r) = w(r)~1 dip(r)/dr is the logarithmic derivative of the wavefunction W(r) and the potential U(r) is a piecewise continuous real 


function of the variable r. By investigating the one-dimensional Schrodinger equation via the exact quantization rule, Qiang and Dong obtained a 
proper quantization rule expressed as 


ie ky(r) dr — cite ko(r) dr = na, (3) 
Roy= zm 0) (4) 


where 1, and rp are two turning points determined by E = U(x). is the number of nodes of #(r) in the region E > U(r) and it is larger by one than 
the number of nodes of wave function w(x). In a similar way, they generalized the quantization rule (3) to three-dimensional Schrodinger equation 
with a spherically symmetric potential. This is realized by taking p(r) = r~1R(r)Y,, (@,@). Substituting this into the Schrodinger equation allows 
us to obtain the following radial Schrodinger equation 


1(1+1)h? 
2mr2 


ZERO) = FEE Ver MIRO, Vere) = VO) + (5) 


from eq.3, the two integrals have the same mathematical form. Accordingly, when applying it to calculate the energy levels, we can calculate its first 
integral with respect to k(r) and then replace energy levels E£,, in the result by the ground state energy Ey) to obtain the second integral. 


QUANTUM MECHANICAL GRAVITATIONAL PLUS HARMONIC OSCILLATOR POTENTIAL 
In this work, we study the one dimensional Schrodinger equation for the QMGHOP model [30], 


V(r) =yr + Voe + smw*r?, (6) 
Eq. (6) can also be written as 

Vir) =V+(v-av)r+ (Gaz + +mw?)) v2, 
Where the turning points r, and r, are determined by E = U(r), 
E=VW+(v-av))r+ (Gaz + +mw?)) rv? (7) 


For simplicity, we can assume the following 


(az +mw*)) = A,(y —aVo) = B,E-Vy =C (8) 
—b-vb?2—4c —b+vb?—4c 
a (9) 
where we have obtained the two turning points 
= Bee _—— Y-QAVo = Ge = E-Vo 
Bete =~ Cay tmnt © a 2 = — Cavs Sie?) 10) 


The momentum k(r) between the two turning points is given by 


k(r) = EC - ar? ~ Br) (11) 


The corresponding non-linear Riccati equation (2) becomes 
d 2 
a PT) = — FF (Co — Ar? = Br) — $(r)? (12) 


Based on the Sturm—Liouville theorem, after taking dp) (r) = —C,r + Cz(C, > 0), substituting this into the above Riccati equation (11) we obtain 
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c= NR _#t 
0 vam 4A 


on considering the first integral in eq.3, we obtained 


Si? ken(r) dr = |” ta — 7) ar 


n° 


_ mv2mA 2 
are (B2 + 4AC,) 


Where we used the following integral formula [31] 
i (r—1)(r2 —1r)dr = 5 (tr —7)? 0<ry <7 


By replacing C,, in (14) by Cp given in (12) and considering (3), we obtain the following result: 


Ga = W2mA (n+ 1) = 


Substituting eq. 8 into eq. 17, we obtain the energy eigen value 


a m 


For the sake of simplicity throughout this study, we let the following 


2m(a2Vo+ 5m?) 


=O, 
m 


(y-aVo)” 


4(a?Vo+ Smw?) id 


E, =(n+5)ho +0 


THERMODYNAMIC PROPERTIES OF THE SYSTEM 


In other to consider the thermodynamic properties of the QMG+HO potential for s-wave state, we calculate first, the vibrational partition function 


Q(B) = Din=0 eBEn, 


B=—_ , Where k is the Boltzmann constant. 


Q1”” (B) = Fxg e Alm +a)>e] 


Bho 


17? (B) =e 2 P? Yee Pan 


F _g(@ Bho Bye 
@,7(B) = e A(Z*) (=) = e-#e (5) 


F -1 
Gg," (8) = eFe (2sinh on 
The vibrational partition function with N-body interaction with no interaction inside is given by 


Qn?(B) = Q1" (B)" = e-NBe (2sinh a 
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Partition FunctionQ(B) 


Figure 1 Vibrational partition functionQ Pid against a for different 6 (0.5, 1, 1.5 and 2) 


Taking natural log to both side of equation (28), we obtain 
InQu”” (B) = Q1" (B)" = —NBo — Nin (2sinh 2) 


The vibrational mean energy U(B) 
The Vibrational mean energy for the QMGHOP model is obtained as 


U(B) = ~ Fr InQy”?(B) 


u(p) =-2 (-Néo — Nin (2sinh 3) 


- ho Bho 
U(f) N(o+ - coth = ) 


(29) 


(30) 


(31) 


(32) 


Mean EnergyU(B) 
Ww 
1 


B 
B 
B 
B 


0. 
1 
1. 
2 


wu 


u 


0.0 0.5 1.0 


Figure 2 Vibrational Mean Energy U against a@ for different 6(0.5,1, 1.5 and 2) 


The vibrational mean free energy F(f): 


F(B) = —FInQu””(B) 
F(g) =- ; [-nBe - Nin (2sinn=)| 


F(B) = [Ne + ain (2sinn)| 
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F(B) = [Ne + NkTIn (2sinh =) (36) 


Mean Free EnergyF(B) 


Figure 3 Vibrational Mean Free Energy F against @ for different 6(0.5, 1,1.5 and 2) 


The Chemical potentialu(N) 


w= = (Ne + NkTIn(2sinh)| (37) 
w= [e + kTIn (2sinh =) (38) 


Chemical potentialu(N) 


Figure 4 Chemical potential w(N) against @ for different 2 (0.5, 1, 1.5 and 2) 


The Vibrational specific heat capacity C,,(B) 


Cy(B) = SP = — Kp? OO (39) 
C,(B) = —kB? a [v (9+ =coth™)| (40) 
C,(B) = Seat (sinn? (m2) (41) 
c,() = 2 (sinn? ()) (42) 
Cy(B) = NkB*y?a? (43) 
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Figure 5 Vibrational Specific heat capacity C, against a@ for different @(0.5,1, 1.5 and 2) 


The Vibrational entropyS (8) 


S(B) = klnQy”(B) — kB é Indy” (B) on 
S(B) =k (-Née — NIn(2sinh fe)) — kB (-N (9+ @coth 3) ae 
S(B) = —Nkin (2sinh™) + kp coth a 


EntropyS(B) 


Figure 5 Vibrational Entropy S against @ for different @(0.5, 1, 1.5 and 2) 


If we let 6 = Vy = 0, then we obtain the energy spectrum of the harmonic oscillator potential. 


En, = (n+5) be 


CONCLUSION 

In this study, we obtained by using the proper quantization rule, the analytical solutions of the Schrédinger equation for the Deformed Harmonic 
Oscillator potential. Furthermore, we study the thermodynamic properties of this model potential. From a high temperature partition function 
obtained, thermodynamic functions, such as vibrational mean energy U, specific heat C, free energy F and vibrational entropy S are derived. It is 
found that in the classical limit, the vibrational entropy S converges to a point for all the diatomic molecules at a different range of B. Also, we have 
conducted the numerical analysis by calculating the bound state energies in terms of the thermodynamic functions. 
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